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^ . Abstract 

' We investigate barotropic perfect fluid cosmologies which admit 

an isotropic singularity. From the General Vorticity Result of Scott, 
it is known that these cosmologies must be irrotational. In this paper 
we prove, using two different methods, that if we make the additional 
assumption that the perfect fluid is shear-free, then the fluid flow must 
CJ , be geodesic. This then implies that the only shear-free, barotropic, 

perfect fluid cosmologies which admit an isotropic singularity are the 



00 

o 



K- FRW models. 

OJQ. 



|xj ■ 1 Introduction 

a' 



In 1985 Goode and Wainwright |[| introduced the concept of the isotropic sin- 
gularity (IS) (see Appendix |A] for the definition) to the field of mathematical 
cosmology in order to clarify what is meant by a "Friedmann-like" singu- 
larity. A question which naturally arises is the following: "precisely what 
cosmologies actually admit an isotropic singularity?" In order to make the 
problem more tractable, the specialisation to perfect fluids has been studied 
by various authors with some success. Goode M has shown that perfect fluid 
cosmologies satisfying the dominant energy condition, with a 7- law equation 
of state, which admit an isotropic singularity, must be irrotational. Scott || 
has extended this result to include all barotropic perfect fluids satisfying the 
dominant energy condition — this is known as the General Vorticity Result 
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(GVR). In this paper we will enlarge the current knowledge of perfect fluid 
cosmologies which admit an isotropic singularity to include the following 
theorem. 

Theorem 1.1 (Zero Acceleration Result) If a space-time (A4,g): 

• is a C 3 solution of the Einstein field equations (EFE), 

• has barotropic perfect fluid source, 
and the unit timelike fluid congruence u is: 

• shear-free, 

• regular at an isotropic singularity (with — 1 < A < 1), 
then the fluid flow is necessarily geodesic. 



We note that the technical condition — 1 < A < 1 implies that the 
barotropic perfect fluid cosmology satisfies the dominant energy condition. 
The dominant energy condition can be used instead of the technical condi- 



tion — 1 < A < 1 in the statement of Theorem |1 . 1| . In the interests of brevity, 
it will henceforth be assumed that when we say "barotropic perfect fluid cos- 
mology" we actually mean "barotropic perfect fluid cosmology satisfying the 
dominant energy condition". 

The General Vorticity Result enables us to say that barotropic perfect 



fluid cosmologies with non-zero vorticity do not admit an IS. Theorem [LI 
which we will refer to as the Zero Acceleration Result (ZAR), implies that 
a shear-free, barotropic perfect fluid cosmology with non-geodesic fluid flow, 
also does not admit an IS. Since the Friedmann- Robert son- Walker (FRW) 
cosmologies are characterised by their globally vanishing vorticity, shear, and 
acceleration, the GVR and ZAR can be combined to produce the fact that 
shear-free, barotropic, perfect fluid cosmologies which are not FRW models, 
do not admit an IS. 

The following conventions are used throughout this paper. 

• Latin letters denote 0,1,2,3. Greek letters denote 1,2,3. 

• * denotes that the entity to which it is attached exists in the unphysical 
space-time (*Ai,*g), not the physical space-time (Ai,g). 

• ; denotes the covariant derivative w.r.t. the physical metric g . 

• *; denotes the covariant derivative w.r.t. the unphysical metric *g . 
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• A dot above an entity in the physical space-time means that we take 
the covariant derivative of the entity w.r.t. the physical metric g in the 
direction of the fluid velocity field, u, i.e., A = A. a u a . Similarly, a dot 
above an entity in the unphysical space-time means that we take the 
covariant derivative of the entity w.r.t. the unphysical metric *g in the 
direction of the unphysical fluid velocity field, *u, i.e., *A = *A* a *u a . 

• A = 0(B) means —K < < K as x — > Xo, for some K £ R + . 

• A = o(B) means — > as x — > Xq. 

• Two functions, A and B, are said to be asymptotically equivalent as 
x — > Xo, written A(x) ~ B(x), if A (a;) = -B(x){i + o(l)} as a; — > x . 



2 Proof of the Zero Acceleration Result by 
elimination 

We proceed to prove the Zero Acceleration Result (ZAR) using two different 
methods. In this section we present a proof by elimination, based on a 
theorem due to Collins and Wainwright 0] . The importance of this method 
of proof, given that we provide a general proof of the ZAR in Section |3], lies 
in the demonstration of the numerous techniques that we now have at our 
disposal for deciding whether or not a specific given space-time actually has 
an IS. We commence by quoting the Collins and Wainwright theorem. 

Theorem 2.1 Any barotropic, irrotational, shear-free, perfect fluid solution 
of the EFE with non-zero expansion and \i + p ^ is either 

1. a FRW model, or 

2. a spherically symmetric Wyman model, or 

3. a plane symmetric Collins- Wainwright model. 
Proof of the ZAR by elimination: 



Scott H has shown (GVR) that a space-time which satisfies the conditions 



of Theorem |T7T] must be irrotational. It is also shown in that if a perfect 
fluid space-time has an isotropic singularity at which the unit timelike fluid 
congruence is regular, then /i + p ^ as the singularity is approached. (More 
precisely, if a perfect fluid space-time has an IS at which the unit timelike 
fluid congruence is regular, then there exists an open neighbourhood U of 
the spacelike hypersurface T = in *M. such that /x + p ^ anywhere in 
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U n M..) Finally, we note that Goode and Wainwright |1| have shown that 
the expansion of the fluid limits to positive infinity as the IS is approached, 
so that the fluid must have non-zero expansion. 



We have now shown that if the conditions of Theorem 1.1 are satisfied, 



then the conditions of Theorem |2.1| are also satisfied. This means that we 
have only three possible models to consider. Both the spherically symmetric 
Wyman models and the plane symmetric Collins- Wainwright models need 
to be eliminated from consideration, since neither has a geodesic fluid flow. 
We will do this by proving below that these models do not, in fact, actually 
admit an IS, and thereby do not satisfy the conditions of Theorem |1.1| . This 
then leaves only the FRW models with their geodesic fluid flows and we are 
done. 



2.1 Spherically symmetric Wyman models 

In local comoving coordinates (t,r,9,(p), the metric of the spherically sym- 
metric Wyman models 0, |5| has the form 



1 2 -dt 2 + dr 2 +r 2 (d0 2 + sin 2 °> 2 



At+B 

where U = U(v ) 7^ 0, v = t + r 2 , a prime denotes differentiation w.r.t. v, and 

U" = U 2 ^ U' 2 = ft/ 3 - \A. (2) 

The quantities A and B are constants satisfying A 2 + B 2 7^ 0, i.e., it is not 
permitted that both A = and B = 0. The energy density, /1, and the 
pressure, p, of the perfect fluid are given by, 

20/7 4 

fi = 3(Av+B) + 12UU', ^ + p=—. (3) 

Let u be the unit flow vector of the fluid, so that g(u, u) = —1. Since we 
are using comoving coordinates, u has the form 

u = u ^-, u° > 0. (4) 
dt v ; 

Thus goou°u° = — 1. Substituting for g 00 we obtain 

U' 2 1 

' ,0\2 



u y = -1 (5) 



U 2 At + B 

* (^°) 2 = ^ (At + B). (6) 
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If a spherically symmetric Wyman model has an IS, then each flow-line 
given by r, 9, <fi constant will encounter the IS at a particular value of t (where 
t could be infinite). That is, t = t(r) is the "equation" for the IS in the 
physical space-time. 

We note, firstly, that the metric will only be Lorentzian when At + B > 0, 
which will henceforth be a requirement. This implies that: 

if A = 0, then B > te (-00,00), (7) 
if A > 0, then t > -f te(-f,oo), (8) 
if A<0, then t < -f te(-oo,-f). (9) 

The expansion scalar # for the fluid is given by 



-3y/At+B if §j > 



Sv/^t+E if ^7 < 0. 



(10) 



It is known |IJ that as we approach an IS along a flow-line, 9 — > +00. 
From the expression for given above it is readily seen, therefore, that metrics 
with A = do not have an IS. Also, for metrics with A > or A < 0, no 
flow-line can encounter an IS at to = t(ro), where to £ since the expansion 
scalar is not infinite there. This leaves only two possibilities: 

1. if A > 0, the IS occurs at t — +00 along every flow-line, or 

2. if A < 0, the IS occurs at t = —00 along every flow-line. 

Since U(v) 7^ 0, then along any given flow-line either U > or U < 0. 
Also U" = U 2 > 0, so that U' is a strictly monotonically increasing function 
of t along each flow-line. 

We firstly consider case ([!]) above. It is known that as we approach 
an IS along a flow-line, fi + p —>■ +00. For the spherically symmetric Wyman 
models, 

u + p= . (11) 

p y 3 U' K 1 

Now suppose that, as t — ► +00: 

a. U' -> -/3 or 0, where /3 G M + . 

This implies that /i + p < as t ^ +00, so that an IS cannot occur 
there. 



5 



U' -> P, where /3 e M + . 
From the equation f/ 2 = 



3 U 



| A, it follows that 



Thus 



20 [§(/? 2 +±A)] 







as i 



-oo, 



;i2) 



(13) 



so that an IS cannot occur there. 



U' 



+oo. 



This implies that U is a strictly monotonically increasing function of t 
along the flow-line as t — > +oo. Indeed, from the constraint equation 
U 2 = |[/ 3 — \A, it may be seen that U — > +oo as t — > +oo. This case 
requires further examination. Using the asymptotic relationship 



U' 



V3 ' 



we obtain the following asymptotic relationships: 

/' 

~ 1 a 



3At + 12 ^§ Ui 

v3 



and/x + p ~ 10^C/2. 



(14) 

(15) 
(16) 



It is known that as we approach an IS with A = — oo, /i = o(/i + p). 
Since, in this case, fi > fi + p &s t — > +oo, an IS with A = — oo cannot 
occur there. 



It is also known || that as we approach an IS with — oo < A < 1, 
fi + p* f(2-A)/x. 
would have to hold: 



jj, -\-p « |(2 — A))tt. In this case, the following asymptotic relationship 



lO^C/l « 2(2- A)At + 8^1(2- X)m 
v 3 v 3 



^ (2-A)At « -W4A-3)[/^ 
v3 



(17) 
(18) 



Now 2 - A > 1 and A > 0. If A < §, then 4A 



3 < 0, and it may be 
seen that the above asymptotic relationship cannot hold. Thus an IS 
with — oo < A < | cannot occur as t — > +oo. 
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If I < A < 1, then 



atxU*, where a x = ^= - G M+ (19) 



1 



2 



* V « [-] f. (20) 
Now [/' « ^ C/i, so that 



3 



U' « fttt, where /5 A = ^ f — 1 " e R+ - < 21 ) 

V3 \ a A/ 

Using L'Hopital's rule yields 

tf«!ftff, (22) 



which is not consistent with the asymptotic relationship ( P0|) given 
above. Thus an IS with | < A < 1 cannot occur as t — > +oo. 

The only remaining possibility is that as t —>■ +oo, we approach an 
IS with A = |. For this case, t = o(U*). It is known ]l| that as we 
approach an IS along a flow-line, 



In this case, 

iabcdr-i _ 256 „4 



;W -> 0. (23) 



C abcd C a6c , = f^f/ 6 , (24) 
(256C/ 6 £/' 2 - 768U 3 U' A + 576C/' 6 ) r 4 



RP9 f> - 



+ (l92r 2 ?7 3 [/ 2 - 288rV 4 - 18A1717' + 216t/[/ 3 



+ 
+ 



u' 2 

(- 864[/[/ 5 - A8AU 4 U' + 72AUU' 3 - 192U 7 U' + 864[/ 4 [/ 3 ) r 2 

^2 

48C/ 8 - 192C/ 5 [/' 2 - A8AU 2 U' 2 + 336t/V 4 + 2AAU 5 + 3^ 2 C/ 2 

[T 72 " 



(25) 
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One can analyse the asymptotic behaviour of R pq R pq using the rela- 
tionship 

_ ^ A 1 

81P 



V3 



(26) 



It is found that the largest terms, namely those of order U 6 , cancel 
one another, as do the next largest terms, namely those of order . 
Therefore, R pq R n = o(ETr). It follows that 



u\ = o[- (27) 



R pq R pg 



Thus the condition that ° Rvq ( ^ bcd — > as t — > +oo does not hold, so 
that an IS with A = | cannot occur there. 

We conclude that for spherically symmetric Wyman models with A > 0, 
no IS occurs at t = +oo. 

Now consider case (|2]), where A < 0. We will examine the possibility that 
an IS occurs at t = — oo along every flow-line. Recall that as we approach 
an IS along a flow-line, fi + p — > +oo. For the spherically symmetric Wyman 
models, 

20 U 4 

< 28 > 

Now suppose that, as t — » — oo: 

a. U' — ► — oo or — j3, where (3 e M + . This implies that fi + p < as 
t — >• — oo, so that an IS cannot occur there. 

b. U' — > (3, where f3 G 1R + . From the equation f/ 2 = |t/ 3 — it follows 
that 

U ^ [U(3 2 + \A)]K (29) 

Thus 

M +p^^Ii^!±Mll e R ast^-oo, (30) 
so that an IS cannot occur there. 

c. U' — > + . This implies that U is a strictly monotonically increasing 
function of t along the flow-line. So we can set U = ^(3A)3 + H(t), 
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where H(t) — > + and H' = U' > as t —>■ — oo. Substituting these 
expressions for U and U' into the equation U 2 = |f/ 3 — jA, we obtain 



(H 



t\2 



l A + ^(3A)§fT + ~(3A)^H 2 + if 3 
8 4 V ' 2 y ' 



Since if = o(l) as £ — > — oo, this implies that, as t — > — oo, 

1 



<£> if' 



(3A)3i? 



v/2 



(3A)3i/2 (since A < 0) 



(31) 



(32) 



(33) 
(34) 

(35) 



Using L'Hopital's rule yields 



2m « - ^=(aA)*t, 



(36) 



which cannot be true since 2if2 — » + but — -^(3A)3t -00 as 
t — > —00. We conclude that the asymptotic behaviour U' — ► + as 
i — > —00 is not consistent with the constraint equation U 2 = ^U 3 — jA, 
which must be satisfied by these models. 

Thus for spherically symmetric Wyman models with A < 0, no IS occurs 
at t — —00. 

This then proves that the spherically symmetric Wyman models do not 
admit an isotropic singularity. 



2.2 Plane symmetric Collins- Wainwright models 

In local comoving coordinates (t,x,y,z), the metric of the plane symmetric 
Collins- Wainwright models M has the form 



IP 



^-dt 2 + dx 2 + e~ 2x {dy 2 + dz 2 ) 
m 2 



(37) 



where U = U(v) 7^ 0, v = t + x, a prime denotes differentiation w.r.t. v, and 

U" + U' = -U 2 . (38) 
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The quantities C and m are positive constants. The energy density, /x, and 
the pressure, p, of the perfect fluid are given by, 



1 2U A 
p = — [3m 2 - 2U 3 - 3(U'+U) 2 ] , (JL + p=—. (39) 

Let u be the unit flow vector of the fluid, so that g(u, u) = —1. Since we 
are using comoving coordinates, u has the form 

u = u° > 0. (40) 

Thus goou°u° — — 1. Substituting for (700 we obtain 

~ -2 772V) 2 = - 1 ( 41 ) 
The expansion scalar 9 for the fluid is given by 

r 3^ if ^>o 

e={ c U TJ (43) 

\ -3f if ^7 < 0. 

The flow-lines of the fluid are given by x = constant, y = constant, and 
z = constant. It is known [1] that as we approach an IS along a flow-line, 
6 — > +00. Since, for the plane symmetric Collins- Wainwright models the 
expansion scalar is either a positive constant or a negative constant along 
each flow-line, it is clear that the condition 9 — ► +00 is never satisfied along 
a flow-line. Thus the plane symmetric Collins- Wainwright models do not 
admit an isotropic singularity. 



3 General proof of the Zero Acceleration Re- 
sult 

The proof by elimination of the ZAR utilises the theorem of Collins and 
Wainwright . Given the existence of the proof by elimination, there should 
also exist a general proof, which we now give. We believe that the general 
proof may also provide direction on how to proceed in proving the FRW 
conjecture (see the discussion in Section [|). 
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General proof of the ZAR: 



In order to prove Theorem we will show that, for a space-time which 
satisfies the conditions of the theorem, the electric and magnetic parts of the 
Weyl tensor must be zero. A barotropic perfect fluid cosmology with zero 
Weyl tensor must be a FRW model, and hence the fluid flow is geodesic. 

The General Vorticity Result of Scott H ensures that a space-time which 
satisfies the conditions of Theorem |LT]must be irrotational — i.e., the vortic- 
ity, u> 2 , of the fluid must satisfy u> 2 = 0. Furthermore, a space-time satisfying 
the conditions of Theorem |1.1| is also shear-free, i.e., a 2 = 0. 

The vorticity, u 2 , and the shear, a 2 , of the fluid are related to the vorticity, 
*u 2 , and the shear, *a 2 , of the unphysical fluid flow, *u, in the unphysical 
space-time (*M,*g) by 



*u 2 = Q 2 u 2 , *a 2 = Q 2 a 2 

*uj 2 = * a 2 = 



,2 n _^ :-2 f\ (44) 



The constraint equation |7], pl30 eqn 4.19] for the magnetic part of the 
unphysical Weyl tensor is given by 

*H ad = 2% {a *u d) - *h at *h ds (*u% c + *a% c ) * V s)fbc *Uf, (45) 

which reduces to 

*H ad = 0. (46) 

The physical and unphysical magnetic parts of the Weyl tensor are related 
by 

^ = *H abj (47) 

and thus the magnetic part of the Weyl tensor for the physical space-time is 
zero as required. 

The Einstein tensor for (Ai,g) can be decomposed relative to the fluid 
flow, u, as follows: 



G ab = Au a u b + Bh ab + S a M fe + S fe M a + S afe , (4£ 



where 



A = G cd u c u d , B = \h cd G cd , 

S a = —h ac R cd u d , S ab = h ac R cd h bd — \h ab R cd h cd 



and 



(49) 

Kb = g a b + U a U b . (50) 
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A perfect fluid cosmology satisfying the conditions of Theorem is 
irrotational, and the unit timelike fluid congruence is regular at an isotropic 
singularity. We have the freedom to choose a cosmic time function T defined 
on and a conformal factor f)(T), such that the unphysical fluid flow *u 
in *A4 is orthogonal to the spacelike hypersurfaces T = constant. Using such 
a cosmic time function, we will employ comoving normal coordinates (T, x a ) 
based on the hypersurface T = in The hypersurface T = in *A4 is 
referred to as the isotropic singularity. 

A perfect fluid satisfies the field equation Yi ab = 0. Thus, in normal 
coordinates ([[I|, pl08], also equation ([55])), the unphysical counterpart *E ab 
is given by the equation 

*S a6 = -2F^*a ab (51) 
=> *Z ab = 0. (52) 

Written in terms of geometrical quantities, the unphysical shear propa- 
gation equation || p29] is given by 

W (*& ac - % {a » c) ) + \*h rs (*u\ a - 2*a 2 - *u 2 ) - ±*E rs 

-*ii r *u s + *a rd *a d s + l*6*a rs + *u r *uj s + *E rs = 0, (53) 

which, in the present situation, reduces to 

- W*«(a*c) + \*Ks*u a ^a - %% + *E rs = 0. (54) 

The (0a) equations are trivial, so setting r = fi,s = v and raising these in- 
dices, and noting that in normal coordinates *h a0 = 0, we obtain the equation 

- W a *h v P\ a *p ) + § - W*u u + *E^ = 0. (55) 

The Bianchi identity dealing with the propagation of the electric part of 
the unphysical Weyl tensor is given |7|, pl31 eqn 4.21d using the Appendix 
pl79 eqn 4.21d] by the equation 

*h a m *h c u E ac + *h a im *r] t)rsd *u r *H a s , d - 2*H q (t *r] m)bpq *u b *u p 
+*h mt *a ab *E ab + *6*E mt - 3*E s {m *a t)s - *E s ( - m *u t)s = 
-i(M+*5)V m - %^ m ^ - l*h ta *h mc *£ {a » c) - \% t a *h m c *t ac 

which, in the present situation, reduces to 

*h a m *h c u E ac + *9*E mt = -%^ m ) - l*h ta *h mc ^ {a , c) + ± {*Z a * a +*u a *Z a ) *h mt . 

(57) 
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Since a perfect fluid satisfies the field equation S a = 0, in normal coordi- 
nates (0, pl08], also equation (|100|) ), the unphysical counterpart *E a is given 
by the equation 

T a = 2F— *ii a , (58) 

and thus the terms on the right hand side of the reduced Bianchi identity 
(equation fl5"T|)) become 

% (t^ m ) = 2F^*u t *ii m (59) 

*E (a%/3) = -2F^*u a *u fi + 2F^\ a *0 ) (60) 

1]% = -2F^*u a *u a + 2F^*u% (61) 

*u a *T, a = 2F^*ii a *ii a . (62) 

For the reduced Bianchi identity, the (0a) equations are trivial, so setting 
m = fi, t = v, the Bianchi identity further simplifies to the equation 

V*h p v *E al3 + WE I1V = F^ - *h av *h^*u {a , p) + \*u a ,/h txv ) . (63) 



Now inserting the shear propagation equation ( |55| ) , and noting that in normal 
coordinates *h a ^ = 5 a ^, we obtain the very simple equation 

*E^ V + *Q*E^ U = —F—*E^ U (64) 

*E tiV + T» 0a *E a " + T v 0a *E a ^ + ( % + ^ j *E» V = 0. (65) 



f n 



In normal coordinates, 



*TM lW n W — F* V~f* (c?f?\ 

1 Oa — 2 9 9-yafi i V a — j <? y 7 a,0 l DD J 

TV = — (67) 



We also know that 



19^ = + (68) 
flV a . (69) 



Thus 



T% Q = ^tfV (70) 
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Inserting this expression for T^ Qa into equation fl65|) yields the equation 

*E^ + — + — )*E flv = 0, (71) 



3F 

from which we readily obtain the equation 

(tt*E^) i0 + ^(STE^) = 0. (72) 

In order to solve this system of partial differential equations, we first need 
to obtain some initial conditions. From Goode and Wainwright [0], pl08] we 
know that 

^ a b \T=o = ^fS a b \ T =o and (73) 

*E a b \ T=0 = ^V|t=o, (74) 

where ^S a b \T=o is the trace-free Ricci tensor of the isotropic singularity (i.e., 
the hypersurface T = in *Ai). Since *E a 6 = 0, it follows that 

%V=o = 0, (75) 

and thus 

(76) 
(77) 
(78) 

Since we are using comoving coordinates, the coordinates x a are constant 
along a flow-line, and thus we can treat the partial differential equations 
as ordinary differential equations, in the variable T, along the flow-line. If 
VL*E^ U was defined on an open neighbourhood of T = 0, we could then use 
a standard existence/uniqueness theorem (see Appendix ||) to show that 
VL*E^ U = 0. However, Q*E f11 ' is only defined on the interval [0,6), where 
6 is some positive constant. Suppose then that Q*E >1U ^ anywhere on 
(0, 6). (If VL*E^ V = at any point in (0, b), then we could apply the standard 
existence/uniqueness theorem to obtain Q*E flu = on [0,6)). Let a G (0,6) 
and set (o) = (lW) a ^ 0. Then 

= - ^ (^£^ by assumption) (79) 
(Q*£^)(T) = {Q*E^) a e~^^ dT . (80) 



*77I b\ 

&a |r=o 


= 


* 77>af> 1 
> Hi |t=0 


= 


(1W)|t=0 


= 
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Now since F and *9 are, respectively, at least C 3 and C 2 and F ^ on an 
open neighbourhood of T = 0, 



|| dr = 0{l) as T^0 + . (81) 
3F 

It follows that e~^& dT -> as T -> + , and since (fi*E^) a ^ 0, 

^4 as T -> + , (82) 

which contradicts the above initial condition (equation fl78"|)). Hence, along 
the flow-line, 

tt*E^ = (83) 

*E' W = since Q > on (0, 6) (84) 

*B a6 = 0. (85) 

The physical and unphysical electric parts of the Weyl tensor are related 

by 

E ab = *E ab (86) 
=> E ab = 0. (87) 

Since E ab = along every flow-line, the electric part of the Weyl tensor for 
the physical space-time is zero, as required. 



4 Discussion 

It is interesting to note that there may exist barotropic perfect fluid cos- 
mological models, with unit timelike fluid congruence which is regular at 
an IS, which do not satisfy the dominant energy condition, yet which are 
irrotational. Such models can be accommodated in an alternative version 



of Theorem |1 . 1|, given as Theorem 4.1. Neither the proof of the ZAR by 



elimination (Section g), nor the general proof of the ZAR (Section f|) use 
the dominant energy condition, except, in so far as it is needed to prove the 
GVR of Scott H . We can therefore remove the dominant energy condition 
assumption from the statement of Theorem [Ll] (i.e., — 1 < A < 1) and in- 
stead replace it with the assumption that the fluid flow is irrotational, and 
the two proofs will proceed exactly as before. 
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Theorem 4.1 (Alternative Zero Acceleration Result) If a space-time 
(M,g): 

• is a C 3 solution of the Einstein field equations (EFE), 

• has barotropic perfect fluid source, 
and the unit timelike fluid congruence u is: 

• irrotational, 

• shear-free, 

• regular at an isotropic singularity, 
then the fluid flow is necessarily geodesic. 

Perfect fluid space-times can be categorised according to their vorticity, 
shear, acceleration, and expansion. Using this categorisation, all currently 
known results about barotropic perfect fluid cosmologies with a unit timelike 
fluid congruence which is regular at an isotropic singularity are summarised 
in Figure [l|. 

The abbreviations used in Figure [I] are: 

GW85 Goode and Wainwright's 1985 paper |l| on Isotropic Singularities 
GVR The General Vorticity Result M due to Scott (note: — 1 < A < 1) 
ZAR The Zero Acceleration Result 
FRWR The FRW Result § 

NO IS The unit timelike fluid congruence is not regular at an isotropic 
singularity 

Figure [I] shows that much is already known about barotropic perfect fluid 
cosmologies in relation to isotropic singularities. Goode and Wainwright [l| 
have shown that the expansion of the fluid must approach infinity at an 
isotropic singularity, and hence, whenever the expansion of the fluid does 
not approach infinity at the initial singularity, there can be no IS. Scott 
has shown that for barotropic perfect fluid cosmologies satisfying the 
dominant energy condition, if the space-time has an IS, then the vorticity 
of the fluid is zero. Hence, whenever the vorticity is non-zero, there can be 
no IS. We have just shown (ZAR) that if the fluid is irrotational, shear-free, 
and not geodesic, then there can be no IS. In a forthcoming paper || we find 
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UNKNOWN 



NO IS 

(ZAR) 



u = 



IS / NO IS 
(FRWR) 



Figure 1: Known results about barotropic perfect fluids 



the precise necessary and sufficient conditions for a Friedmann-Robertson- 
Walker model to admit an IS. The only remaining unknown cases, therefore, 
are for irrotational perfect fluids with non-zero shear. 

For the case of a barotropic, irrotational, perfect fluid cosmology with 
non-zero shear, there are examples of cosmologies with, and without, geodesic 
fluid flow, both with, and without, isotropic singularities. For example, the 
Kantowski-Sachs models || are irrotational, have non-zero shear and geodesic 
fluid flow and have an IS The cosmologies given by Mars QTOJ are irro- 
tational, have non-zero shear, and are not geodesic, some of which have an 
IS. A review of known examples of cosmologies which admit an IS is given 
in 
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An idea which has often been linked with the concept of an isotropic singu- 
larity 



1121 , [13], [14], [15| is the Weyl curvature hypothesis |fL6| . Non-rigorously 
stated, the Weyl curvature hypothesis says that "the natural thermodynamic 
boundary condition for the universe is that the Weyl tensor should vanish 
at any initial singularity" . The examples given in the previous paragraph do 
not satisfy the Weyl curvature hypothesis which leads us to what is known 
as the FRW conjecture §|,|0. 
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Conjecture 4.1 (FRW) If a space-time (A4,g) is: 

1. a C 3 solution of the Einstein field equations with a barotropic perfect 
fluid source, and 

2. the unit timelike fluid congruence u is regular at an isotropic singularity 
(with — 1 < A < 1), and 

3. the Weyl curvature hypothesis holds, 

then the space-time is necessarily a Friedmann-Robertson-Walker model. 
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Appendices 

A Definition of an Isotropic Singularity 

In 1985 Goode and Wainwright introduced the concept of the isotropic 
singularity (IS). Scott Jl7|, [3| has amended their original definition to re- 
move some inherent redundancy, and this amended definition of an isotropic 
singularity is given in Definition [A.l| . 

Definition A.l (Isotropic singularity) A space-time (A4,g) is said to 
admit an isotropic singularity if there exists a space-time (*Ai,*g), a smooth 
cosmic time function T defined on *Ai, and a conformal factor Q(T) which 
satisfy 

1. M. is the open submanifold T > 0, 

2. g = Q 2 (T)*g on M., with *g regular (at least C 3 and non- degenerate) 
on an open neighbourhood ofT = 0, 

3. Q(0) = and 3 b > such that Q e C°[0, b] n C 3 (0, b] and Q(0, b] > 0, 

4- A = lim r _ >0 + L(T) exists, A ^ 1, where L = i^-(§r) 2 and a prime 
denotes differentiation with respect to T. 
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Definition A. 2 (Unphysical flow) With any unit timelike congruence u 
in M. we can associate a unit timelike congruence *u in *M. such that 

V = Vtu a in M. (88) 

(a) If we can choose *u to be regular (at least C 3 ) on an open neighbourhood 
ofT = in *M. , we say that u is regular at the isotropic singularity, and 

(b) if, in addition, *u is orthogonal to T = 0, we say that u is orthogonal to 
the isotropic singularity. 



B Existence/Uniqueness Theorem 



The following existence/uniqueness theorem can be found in [|18 



Theorem B.l (Existence/uniqueness) If the functions p and g are con- 
tinuous on an open interval a < x < (3 containing the point x = xq, then 
there exists a unique function y = <j)(x) that satisfies the differential equation 

y' + p(x)y = g{x) (89) 

for a < x < (3, and that also satisfies the initial condition 

y(xo) = yo (90) 

where y is an arbitrary prescribed initial value. 



C Comoving Normal Coordinates 

We set up comoving normal coordinates (T, x a ) in (*A4,*g), based on the 
hypersurface T = (see Section 3). Thus 



*g 0a = o, (91) 



which implies that 

Also 
Then, 
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TV = (92) 



(93) 

*C = o 

*h\ = *h a = (94) 
*hj = 6 J 



*u = *u° = 

u a p u y p 

*E aQ = (96) 
*0 aO = (97) 

Since the physical space-time has a perfect fluid source then Q , 

*E ab = -2F^*a ab} (99) 
Q' 

T, a = 2F—*u a . (100) 

u L 
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